Decoherence of a two-qubit system with a variable bath coupling operator 
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We examine the decoherence of an asymmetric two-qubit system that is coupled via a tunable 
interaction term to a common bath or two individual baths of harmonic oscillators. The dissipative 
dynamics are evaluated using the Bloch-Redheld formalism. It is shown that the behaviour of the 
decoherence effects is affected mostly by different symmetries between the qubit operator which 
is coupled to the environment and temperature, whereas the differences between the two bath 
configurations are very small. Moreover, it is elaborated that small imperfections of the qubit 
parameters do not lead to a drastic enhancement of the decoherence rates. 
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I. INTRODUCTION 



Quantum computation provides a substantial speedup 
for several important computational tasks 0, S S 13 • A 
general quantum bit (qubit) consists of a two level quan- 
tum system with a controllable Hamiltonian of sufficient 
generality to implement a universal set of quantum logic 
gates Q. From such a set, an arbitrary quantum algo- 
rithm can be implemented to any desired accuracy lim- 
ited only by decoherence. A universal two qubit system 
requires just single-qubit rotations and one additional en- 
tangling two qubit gate. One important example is the 
controllcd-NOT (cnot) gate which switches the state of 
the second qubit depending on the state of the first qubit 
and is given in spinor representation of a« <g> by 



U, 



CNOT 
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(1) 



up to a global phase. 

Superconducting Josephson charge and persistent cur- 
rent (flux) qubits have been shown to possess the neces- 
sary properties £| to act as a quantum bit. They have 
been manipulated coherently and coherence times in the 
fis range have been demonstrated 0, IS EJ EU EH E3 
with a corresponding quality factor of quantum coherence 
of up to Q v « 10 4 [l(j- In a two qubit system, where the 
coupling was achieved using a shared Josephson junction, 
coherent Rabi oscillations between states of a coupled 
qubit system were observed 0, Q] and in a two charge 
qubit system a conditional gate operation was performed 
|l2j . All of these experiments suffer from material im- 
perfections which lead to non- ideal time evolutions of the 
quantum states due to a parameter spread in the char- 
acteristic energies of the system Hamiltonian. On the 
other hand, for high symmetry of the qubit parameters, 
the qubit coherence can be strongly protected. This ex- 
tends from the protection of the singlet in a symmetric 



qubit setup to t he g eneral concept of decoherence- 
free subspaces (DFS) [16tll7| . General considerations on 
the stability of such DFS can be found in Ref. ^S In 
this paper, the dependence of the decoherence rates and 
gate quality factors on the parameter spread of the qubits 
will be elaborated theoretically. In perspective, this is of 
crucial importance for connecting the experimental sta- 
tus and prospects to these central concepts in quantum 
information science: which degree of parameter unifor- 
mity do experiments have to achieve for symmetry-based 
protection schemes to work — do these schemes have to 
be extended in order to accomodate experimental restric- 
tions? 

Moreover, variable bath couplings to the decohe ring 
environment have also been identified, e.g., in Ref. |l9j| . 
The decoherence properties of different qubit operators 
that are coupled to the environmental bath will be inves- 
tigated in this work. 

The Hamiltonian of a typical pseudo-spin system can 
be expressed in terms of the Pauli matricies as 



H q = 4(e<r 



(2) 



In a two qubit system, an additional interaction term is 
required to implement the universal two qubit g ate. In 
superconducting implementations jo, B S US* UJS l2u UM 
this coupling term is typically proportional to cri^ <ri 2 \ 
Here, the superscripts are the qubit indices. In partic- 
ular, inductively coupled flux qubits [TM lz| and capac- 
itively coupled charge qubits |23j are coupled this way. 
Thus, the two qubit Hamiltonian is 



H 



2q 
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(3) 
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In the sing let/triplet basis, (1,0,0,0) T = |H), 
(0,1,0,0) T - (|U) + |it»/V2, (0,0,1,0) T = |U), 
(0,0,0, 1) T = (|t|) - |IT))/\/2 that exhibits the symme- 
try properties of the coupling most clearly, this Hamilto- 
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nian takes the following explicit matrix form 

i / 
H 2q = -- 



V 




(4) 



with e = + e (2 \ Ae = e« - e< 2 > and V = faCO + 
?7 (2) )/\/2, Ar/ = (j/ 1 ) -7/ 2 ))/\/2. Using this Hamiltonian 
the CNOT gate can be implemented through a sequence 
of elementary quantum gates [TEl |24| 

t/cNOT = 4 2) cxp(- i ^«)exp(-^Jaf) 

{-^^^)exp(-^^))t/f,(5) 



x exp 1 



(2) 

where U# denotes the Hadamard gate operation per- 
formed on the second qubit. It involves one two-qubit op- 
eration at step three only. For our numerical calculations 
we applied the characteristic energies that were used in 
|l5j | as a viable example for superconducting solid-state 
flux or charge qubits. Following this approach also gate 
sequences optimized with respect to decoherence have 
been studied |25j . 

Disregarding the Hadamard gates, the gate operation 
Eqn. (J5J forms a controlled-phase (cphase) gate 



t^CPHASE 
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(6) 



With (f) = 7T. 

In experimental realizations of this model, additional ef- 
fects always impair the capability of the system to op- 
erate as a qubit. In condensed matter implementations, 
the most pronounced is the coupling to environmental de- 
grees of freedom. This leads to relaxation, i.e., classical 
thcrmalization of the states as well as, on a much shorter 
timescale, to dephasing. Decoherence causes the system 
to act like a classical ensemble eliminating all potential 
computational benefits of quantum algorithms. For a 
wide range {e.g. [TEL I23L |2(J) of solid state implemen- 
tations the dominant decoherence effects caused by cou- 
pling to linear environments such as electric circuits obey 
Gaussian statistics and can be effectively modeled with 
a bath of harmonic oscillators. It is assumed here that 
there is only one decoherence source in the dominating 
order of magnitude in the coupling parameter and pos- 
sible weaker noise sources are ignored. To model this 
source each qubit is either coupled to an individual or 
to a common bath of harmonic oscillators. The system 
Hamiltonian then takes the form 

H£ b = H 2g + \ (*« + X^) + Hg> + Hg> 

(7) 



or 



H 2qb - H 2<? I" j 



\ (oM + sw) x + n B , (8) 



where a s is the spin-representation of the qubit opera- 
tor talking to the environment that depends on the spe- 
cific implementation of the qubit. For the special case of 
superconducting flux qubits which only experience flux 
noise and superconducting charge qubits which are only 
subject to charge noise this would correspond to a s = <J Z . 
Here, X is the collective coordinate of the harmonic os- 
cillator bath and the superscript distinguishes between 
the single bath and the two bath case. The general form 
is 



<t s = (c • a) = V2(c x a x + c y a y ) + c z a z , 



(9) 



where the factor V2 in front of c x and c y was chosen 
for convenience in the singlet/triplet basis in which the 
qubit-bath interaction becomes 



H; 



c z X c-X -c-AX 

c+X c^X c z AX 

c+X -c z X c+AX 

-c+AX c z AX c_AI 



(10) 

with c± = c x ± icy, X = + X^\ Here, AX = 

i-(i) _ jy(2) for t he case of two baths and AX = for 
one common bath. 

In the following we will, without loss of generality, 
characterize the results by the angle 9 between the a x 
and a z component of the coupling 



a s X = (a x sin 8 + a z cos6*)X. 



(11) 



This is completely analogous to the bath coupling that is 
encountered in proposed experimental qubit realizations, 
e.g., for charge qubits [l^. The bath coupling angle 9 is 
defined for 6 £ [0,tt/2]. 

Following the lines of 0] and j2|| the Bloch-Redfield 
formalism is applied to calculate the effects of decoher- 
ence. The Bloch-Redfield equations and decoherence 
rates are given analytically. However, in comparsion to a 
fully analytic evaluation of the dynamics of the two qubit 
system |27|. with this method the time evolution of the 
reduced density matrix can also be determined numeri- 
cally for a wide range of system Hamiltonians. 

The environment, i.e., the bath of harmonic oscillators, 
is characterized by its spectral density. The strength of 
the dissipative effects is given by the dimensionless pa- 
rameter a. The bath spectral function is assumed to be 
linear in frequency up to a cutoff frequency uj c . Thus, 
J{uj) — ahbj j '(1 + (w/w c ) 2 ), i.e., we employ an Ohmic 
spectrum with a Drude cutoff. The cutoff frequency is 
chosen two orders of magnitude above the largest fre- 
quency which is typical for a flux qubit system, oj c = 10 13 
Hz 0. 

We choose a rather large coupling strength to the envi- 
ronment of a — 10~ 3 , which is still in the weak coupling 
regime, to be able to observe pronounced decoherence 
effects. 

The Bloch-Redfield equations describe the evolution of 
the density matrix in the eigenbasis of the unperturbed 
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Hamiltonian 28, 2 



here V denotes the principal value. Likewise, 



- 'y^^Rnmkipkt, (12) 
kl 

where the Redfield tensor R nm ke is given by 

J^nmkl — vim L nrrk ~<~ "nk L trrm L Imnk L Iranki 
r r 

(13) 

and the rates T are given by the Golden Rule expressions 



T hLk = K~ 2 I dte-^iHj^Hj.nkiO))? (14) 
Jo 



1 



x (coth(/3?iw£ m /2) + 1) 

A 2 M 2 ~(w im ) + A 1 Aff (w im ) 



(20) 



The limit of ujg m tending towards zero can be evaluated 
separately 

r (+) =r (-) = JL( & W Q.0-) _j_ 
imnk imnk A/31% s ; n & s,im s,nk 1 

+ + ( 21 ) 



^ -ft" 2 



dt e -^" fc *> 



Hi 4m (0)H Iink (t))p, (15) 



where Hi t t m (t) is the matrix element of the bath/system 
coupling part of the Hamiltonian in the interaction pic- 
ture and in the eigenbasis of the decoupled Hamiltonian. 
Here, (3 indicates averaging over the degrees of freedom 
of the thermal bath. In the following (3 = l/fc^T, where 
T is the temperature of the bath. Evaluating this, we 
find according to 0] for one common bath the rates 



r£L = ^AJK fe ) (cothG8fiw nfc /2) - 1) + ^ 



xV / duj 



(coth(f3hLu/2)u> n k — ui) . 



(16) 



with A = K tmnk = ^} m (T^l k + af} m afl k +af} m (T^l h ■ 
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tin 



(17) 



For two distinct baths 



1 



[A'JltWni) + A 2 J 2 (w„ fe )] 



" imnk 
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A^M+K^+A^+Kfe) 



,(18) 



with A 1 — A 1 -rr (1) rr (1) A 2 -A 2 - n- (2) rr (2) 

witn A - A lmnk - <J s £m <J s nk , A - A imnk - cr s lrn a s nk 
and 



oo 

M^(Q)=V J dcj 



j 2 -O 2 



(coth(/3^/2)fi=Fw), (19) 



for one bath, and 
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s,zm s,nk 



a2<7 



(2) ~(2) 



s.im 



^s,nk ) 

^ . ( 22 ) 

for two baths. All calculations we re p erformed in the 
same parameter regime as in Ref. [15J, thus renormal- 
ization effects of the frequencies are weak and will be 
neglected. 

The ability of a realistic circuit, or in our case a more 
realistic model of a circuit, to operate as a quantum 
gate is characterized by the four gate quality factors in- 
troduced in jsjj. Those are the fidelity !F, purity V, 
quantum degree Q and entanglement capability C. The 
quantum degree and entanglement capability character- 
ize entangling operations. They are unique to multi-qubit 
gates. We will collectively refer to these as non-local gate 
quality factors (GQFs) as opposed to fidelity and purity, 
which are both well defined for an arbitrary number of 
qubits, in particular also for a single qubit, and will be 
referred to as the local gate quality factors. 

The fidelity can be evaluated |2J| as follows 



r = {*in\utp aa tU\* t 



(23) 



The overline indicates the average over a discrete set of 
unentangled input states |^i n ) that can serve as a basis 
for all possible input density matrices. The propagator U 
is the ideal unitary evolution of the desired gate, and p out 
is the density matrix after applying the realistic gate to 
| 'fin). Thus a perfect gate reaches a fidelity of unity and 
the deviation from unity indicates the deviation from the 
ideal process. The purity V is indicative of the decoher- 
encc effects, 



V = tr(p 2 out ). 



(24) 



Again, the overline indicates the input state average. A 
pure output state leads to V equal to one, whereas as 
the state becomes increasingly mixed, the square of the 
weight of the contributions no longer sums up to unity 
and goes down to a minimum of one divided by the di- 
mension of the Hilbert space of the system, 1 /4 in our 
case. 
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CNOT 


local GQFs 


nonlocal GQFs 


preferred case 


T<T S 


close to (7z 


close to a z 


1 bath 


T>T S 

CPHASE 


close to a z 
local GQFs 


at (Ta; 

nonlocal GQFs 


1 bath 
preferred case 


T<T S 


at o" z 


at <r z 


a 




at cr z 


at CTjc 


2 baths 



a There is no clear tendency observed in this case, see Fig. □ Close 
to pure u z coupling two baths are preferred, and close to (7 x OI16 
bath. 

TABLE I: Maxima of the gate quality factors for the CNOT 
and the CPHASE gate operation. Here T indicates the tem- 
perature and Ts = Es/ks is the characteristic temperature 
scale, which corresponds to the qubit energy scale during the 
gate operation. Both the preferred bath configuration and 
qubit operator coupling to the bath are given. 

Whereas the preceding two factors can be defined for any 
number of qubits, the following two are particular to the 
higher dimensional case 

Q= max (* me |pout |*mc) (25) 

Pout,|*mc> 

Here the p ou t are the density operators after the gate op- 
eration relating to unentangled input states, whereas the 
| "Pmc) are the maximally entangled states, also known as 
Bell states. Therefore, this measures the ability of the 
gate to create quantum entanglement. 
Finally, C is the smallest eigenvalue of the density ma- 
trix resulting from transposing the density matrix of one 
qubit. As shown in :31], the non-negativity of this small- 
est eigenvalue is a necessary condition for the separability 
of the density matrix into two unentangled systems. Af- 
ter separation, the partially transposed density matrix is 
a valid density matrix as well, with non-negative eigen- 
values. The negativity of the smallest eigenvalue thus 
indicates that the states are not separable and thus non- 
local. It approaches -0.5 for the ideal CNOT gate, the 
dynamics of entanglement in a two-qubit system have 
been studied in Ref. . The entanglement capability 
is closely related to the negativity En of a state |3j| which 
is a non-entropic entanglement monotone |34| . 

II. COMBINATION OF & x AND <r z ERRORS 

Now, the spin-boson model with a variable coupling 
operator to the harmonic oscillator baths is studied 
in more detail. 

We start with the CPHASE gate which is entangling and 
forms the core part of the CNOT operation. The quantum 
degree for the CPHASE gate is always smaller than the 
ideal value because the cphase gate cannot create entan- 
gled Bell states. Thus, we did not consider the quantum 



degree for the CPHASE gate. The different error coupling 
configurations achieve the best gate quality factors for 
different coupling operators to the environmental baths. 
The scenarios are summarized in Table 1. 

Two qualitatively different temperature regions are 
found, separated by a smooth crossover. For low tem- 
peratures T < T s , where E s /h = (k B /h)T s = 1 GHz is 
the characteristic energy scale, which corresponds to the 
typical qubit energy scale during the quantum gate oper- 
ation, spontaneous emission processes dominate. When 
Ts is approached, thermal effects become important and 
for T ss Ts temperature is the dominating energy scale 
as will be discussed in more detail below. 

The CPHASE gate, for pure a z coupling, is protected by 
symmetry because the gate operation and the coupling 
to the bath commute. As was shown previously, all dis- 
turbances vanish here in the limit of low temperatures. 
In this case, spontaneous emission processes are the dom- 
inating decoherence mechanism because absorption and 
excitation processes are effectively suppressed due to the 
bath coupling (i.e. matrix elements for theses processes 
are restricted due to symmetries of the bath coupling 
operators) and the temperature. In this case, the low 
temperature regime can be referred to as the emission 
limited regime. 

The additional a x operation in the CNOT gate Eqn. (JSJ) 
during the single qubit Hadamard operations leads to 
nonvanishing decoherence rates even in the low temper- 
ature limit. The reason is again the competition of pure 
dephasing with emission and absorption processes which 
show a different dependence on the coupling angle 9. 

However, the Hadamard part of the CNOT operation 
is short compared to the overall length of the gate op- 
eration (jSJ). Thus, it is found that for low temperatures 
the best values for the GQFs are obtained very close to 
pure a z coupling as depicted in Fig.^ This implies that 
the overall decoherence effects are smallest if the bath 
coupling angle 6 throughout the gate is adjusted to the 
distribution of gate operations, which are characterized 
by different directions on the Bloch sphere. 

For the CNOT gate better results are observed for the 
single bath case throughout the low temperature regime, 
see Fig.n For the single bath configuration, close to pure 
a x coupling to the bath the difference becomes quite sig- 
nificant and for the non local GQFs actually approaches 
a factor of two, but reduces again as pure a x coupling is 
reached. The CPHASE gate (cf. Fig.^l prefers a two bath 
configuration unless the coupling is o x dominated. Again 
the non-local GQFs are affected most. 

For the pure a z case (9 = or correspondingly c x — 
c y = and c z = 1) little difference between the two bath 
and the single bath behavior is found in the CNOT case 
and none at all for the CPHASE. It is observed that the 
single bath configuration is certainly preferred as soon 
as there is a significant a x contribution in the gate op- 
eration. This means that the additional protection from 
the one dimensional decoherence free subspace |l5lll6j ] in- 
volved is mainly beneficial if the commutator of the qubit 
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FIG. 1: Dependence of the gate quality factors on the bath coupling angle 9 defined in 1111 for the CNOT and CPHASE operation 
at T m <C Ts. Here, the behaviour of the gate quality factors for both the single bath and two bath case is shown. The 
characteristic energy scale for the gate operation is Es/h = 1 GHz |15jl . The lines are provided as guides to the eye. 



operator, which couples to the bath and the qubit Hamil- 
tonian (the Hamiltonian that is needed to perform the in- 
dividual parts of the gate operation) has appropriate ma- 
trix elements, i.e., if there is a significant non-commuting 
part in the bath coupling and the gate operations. How- 
ever, in a <7 2 dominated case the individual coupling is 
preferred as it does not induce any additional indirect 
couplings between the qubits. It is natural that the two 
qubit GQFs should notice this more strongly than the 
single qubit GQFs for which the differences never become 
more than about one fifth of the individual deviations. 
For the high temperature regime, drastically different be- 
havior in the non-local GQFs is found, see Fig. [S] Both 
of them now achieve their best values at a pure a x cou- 
pling for both gates, the local GQFs achieve their maxi- 
mum at a pure a z value. The protection that the CPHASE 
gate enjoyed in the low temperature regime breaks down 
here. The high temperature case is essentially scale free, 
i.e., high temperatures symmetrize the system. In this 
case the system eigenbasis is given by the qubit opera- 
tor which couples to the bath. This can be nicely shown 
when considering the single qubit dephasing rates within 
the spin-boson model |2l|. 



1 

T~2 



2E 2 



5(0) 



2E 2 



S(E). 



(26) 



This becomes w 2irakBT '/% for T 3> E and does not 
depend on the ratio A/e. 

Thermalization is determined by the off-diagonal bath 
couplings in the basis of the corresponding system Hamil- 
tonian which is required for a certain gate operation. It 
will be strongly dominated by the non-commuting contri- 
butions, i.e., the a x bath coupling for the CPHASE part 
of the CNOT gate. The single-qubit Hadamard part of 
the CNOT gate will be additionally also affected by the 
a z bath coupling. The two-qubit operation (CPHASE), 
or in other words the non-local part, of the CNOT gate 
is of the dg ® &z -type and the single-qubit Hadamard 
gates contain both a x and a z contributions. Thus, dur- 
ing the gate operation the thermalization is dominated 
strongly by the a x part of the bath coupling for the 
non-local and by the a z part for the local GQFs, imply- 
ing that for thermal fluctuations the <7 T -typc couplings 
are more important in inducing inter-qubit transitions, 
while a z primarily affects the single qubit gate quanti- 
fiers. What implementation to choose for a gate here, 
becomes a question of what gate quantifiers are desired 
to be optimized. The differences between the one bath 
and two bath scenario are now small. 

For pure a z coupling of the qubit to the bath, a pecu- 
liar effect is observed. In this case, the minimal eigen- 
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FIG. 2: Dependence of the gate quality factors on the bath coupling angle 6 denned in Eqn. for the CNOT operation at 

large temperatures T = 2Ts and T = 200Ts . The characteristic energy scale for the gate operation is Es /h = 1 GHz iHf . The 
lines are provided as guides to the eye. 



value of the partially transposed density matrix that 
is the entanglement capability stays negative even for 
T > T s (Fig. [21). The negativity of the eigenvalue of the 
partially transposed density matrix is not just a neces- 
sary but also a sufficient criterium for the nonseparabil- 
ity of the system in our case [3l]]. Thus, no matter what 
the temperature or strength of the dissipative effects in 
our system during the CNOT gate operation, entangle- 
ment will never be eliminated completely. This can be 
explained by fast thermalization into a protected entan- 
gled state. Furthermore, this effect carries over well into 
the regime where both a x and a z noise are present. 

Overall, the temperature (and the coupling strength) 
has the largest influence on the GQFs. At T -C T s 
(T w 10 -3 K), we observe deviations of the GQFs from 
the ideal value which are less than 10 -3 . At T T s 
(T w 1 K) the deviations are 10 _1 and quickly increas- 
ing further at larger temperatures. The different coupling 
operators to the bath are the next strongest effect. Ro- 
tating the coupling operator from a z to a x causes, in the 
worst cases, three to four times stronger deviations from 
the ideal value than the a z noise. Finally, the change 
due to different types of bath couplings is generally small 



compared to the differences between a x and a z type cou- 
pling. This also suggests that we do not need to worry 
about noise sources with a weaker coupling strength, even 
if they couple through a less favorable coupling operator. 

As an intermediate conclusion, it is found that for 
the decoherence dominated regime the CPHASE opera- 
tion reached the optimum value of the GQFs for a pure 
Decoupling to the bath. In the case of the CNOT oper- 
ation, the minimum is located slightly shifted to the a x 
component because of the mixture of a x parts during the 
Hadamard operations, compare with Fig.^ 

For the CNOT operation, the optimum values of the 
four gate quantifiers are encountered at different bath 
couplings, which are characterized by the mixing angle 
Eqn. especially for large temperatures. 

Thus, the differences between the case of one common 
bath and two baths are much less important than the 
symmetries between the gate operation and the bath op- 
erators. In particular, the difference between the case of 
one or two baths disappears for pure a z coupling to the 
bath. Here, decoherence due to flux noise or charge noise 
in coupled superconducting flux or charge qubits was ex- 
plored. Decoherence due to 1/f noise, caused by back- 
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FIG. 3: Dependence of the decoherence rates on the qubit 
asymmetry. Here, we set K — 0, e, = 0, Ai = Es, and 
vary A2. The single and two bath cases behave identically, 
thus only the single bath case is shown. The strength of the 
decoherence effects is set to a = 10 -3 and T ~ 0.5Tg. We 
set the bath coupling angles 1271 to 9\ — and 62 = 0. The 
decoherence rates are scaled by avs, with vs = Es/h. The 
lines are provided as guides to the eye. 



ground charges or bistable fluctuators, was not treated. 



III. NON-IDENTICAL QUBITS 

Now, we do not restrict the analysis to the case of a 
uniform error coupling (|llf> anymore. In general, both 
qubits can couple to the baths differently 

<xW = df> sin Oi + of> cos 6 l: (27) 

where i = 1,2 denotes one of the qubits. For the nu- 
merical calculations, the qubits are set to the degener- 
acy point K = e, ; = and Ai is set to Es, i.e., ef- 
fectively this model describes a system of two uncou- 
pled qubits. Here, A2 and thus also the asymmetry 
G = (Ai - A 2 )/(Ai + A 2 ) is varied. 

Experimentally the spread of the qubit parameters due 
to fabricational imprecision is very important both be- 
cause quantum algorithms (without further modification) 
require a certain level of precision and the decoherence 
effects in the system of qubits have to be sufficiently small 

m. 

Therefore, it is of central importance to investigate also 
the effects of the parameter spread in non-identical qubits 
on the behaviour of the decoherence rates. Supercon- 
ducting qubits are preferably operated at the degeneracy 
point where decoherence effects are suppressed for super- 
conducting charge and flux qubits. However, the tunnel 
matrix elements for superconducting qubits can differ sig- 
nificantly, on the order of several percent [T2I fl4| . Thus, 
the dependence of the decoherence rates, i.e., dephas- 



FIG. 4: Dependence of the decoherence rates on the qubit 
asymmetry at T w 0.5Ts. Here, the case of one common 
bath is investigated. The tunnel matrix element of the second 
qubit and the inter-qubit coupling are set to A2 = K = Es 
and Ai is varied. For comparison with experiments, large 
asymmetry in the tunnel matrix elements of the individual 
qubits is investigated. The bath coupling angles are set to 
9i = and 62 = 0. The strength of the dissipative effects is 
a = 1CT 3 . The lines are provided as guides to the eye. 

ing rates and the relaxation rate, close to the degeneracy 
point on the qubit asymmetry is an important property. 

Figure |31 depicts the dependence of the decoherence 
rates on the qubit asymmetry G when the individual 
qubits are operated close to the degeneracy point. Tem- 
peratures typical for experimental situations are chosen 
for this analysis, i.e., T w 0.5Tg. 

We observe that for pure a x coupling to the bath 
(81 = tt/2 and 6> 2 = tt/2), the asymmetry of the qubits 
is irrelevant because the coupling to the bath and the 
system Hamiltonian commute (the indices of qubit one 
and qubit two could be exchanged without changing the 
system). For mixed bath coupling of the a x type for one 
qubit (61 = 7r/2) and the a z (# 2 = 0) type for the other 
qubit, still the decoherence rates do not vary for differ- 
ent asymmetry. The reason for this behaviour is that the 
a x bath coupling of the second qubit always commutes 
with the qubit Hamiltonian, i.e., only flipless dephasing 
processes contribute to the decoherence rates of the sec- 
ond qubit. When we vary the asymmetry, essentially A 2 
is varied (Ai = Es is kept constant), which leads to a 
different contribution of the second qubit to the overall 
decoherence. However, these corrections are small com- 
pared to the full decoherence rates (i.e., not only flipless 
dephasing processes) that contribute in the case of qubit 
one where [Hsb, H sys ] ^ and Ai stays constant. 

Finally, Fig-El shows the case of fully perpendicular- 
system Hamiltonian and bath coupling. Here, the deco- 
herence rates increase steeply for increasing asymmetry. 
Note here that due to the definition of G and Ai, the 
two cases G = -0.8 -> A 2 = 9E S and G = 0.8 -> A 2 = 
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FIG. 5: Temperature dependence of selected decoherence rates for K = 0, ej = 0, Ai = Eg, A2 = 0.9Ai, for the single bath 
case. Here, temperature and the bath coupling angles are varied. The strength of the dissipative effects is set to a — 10~ 3 . 
The decoherence rates are scaled by avs, where us = Es/h. The lines are provided as guides to the eye. 



LIE'S are vastly different. 

From the dependence of the decoherence rates on the 
asymmetry G of the two qubits at the degeneracy point 
it is possible to estimate the maximum tolerable asym- 
metry for a given constraint on the deviation of the de- 
coherence rates from their value for perfectly symmetric 
qubits. We define the deviation of the decoherence, i.e., 
the relaxation or dephasing rates from their values at the 
degeneracy point as 

r Rm } = 1 _ r fl,y, 3 ( A i ^ A 2 ,/f = e. t = 0) 
dov r Jl)W (Ai = A 2l A' = e i = 0)" 1 ' 

We find in the case of a bath coupling parallel to the 
Hamiltonian that for r^ ev to be smaller than 1%, it is 
required that O.5A2 < Ai < I.5A2, i.e., the parameter 
spread of the two qubits could be remarkably large (~ 
50%) without considerably affecting the relaxation rates. 
However, detailed analysis shows that for the r?" , i.e. 
the deviations of the dephasing rates, the increase hap- 
pens much earlier. E.g., in the case of already for a 



deviation between Ai and A2 of less than approximately 
10%. Moreover, there is a large spread among the de- 
phasing rates, which are sensitive to the qubit asymme- 
tries. Note that both the single bath and the two bath 
case behave identically for the relaxation rate. Differ- 
ences between the two cases only occur for the dephasing 
rates. The angles of the bath coupling where the min- 
imum dephasing rates are encountered are different for 
the different dephasing rates. 

Typical experimental values for charge 01 an d flux 
qubit designs indicate that the parameter spread in 
the tunnel matrix amplitudes can be quite large, in the 
case of the charge qubit it is a factor of A2/A1 w 0.91 
and for the flux qubit A2/A1 4.22. This difference of 
asymmetries is due to the fact, that fabrication parame- 
ters such as E c and Ej enter the exponent of the tunnel 
splitting in the flux qubit case [(|. These experimental 
values emphasize that it is important to study the evo- 
lution of the decoherence effects for non-identical qubit 
parameters. Moreover, important information about the 
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noise sources coupling to the qubit can be identified. 
From comparison of the decoherence rates for different 
qubit samples, which possess different asymmetries be- 
tween the tunnel amplitudes, it is thus possible to iden- 
tify the predominant bath coupling angle. In most qubit 
designs the bath coupling angle is then uniquely related 
to a certain noise source, e.g., flux noise in the case of 
flux qubits [3(| . 

Figure 0] depicts the experimentally important 0] be- 
haviour of the decoherence rates when A2 and K are 
fixed to Eg and A x > A 2 is changed. Here, different 
from the definition in Eqn. (|28[) . all rates are normalized 
by r^, Vij (Ai = A 2 = K, €i = 0), which will be labelled 
by Tdcv.K- In this case the two qubits are permanently 
coupled, and embedded into one common environmental 
bath. The decoherence rates begin to increase linearly 
when Ai is larger than A 2 . 

Figure [SJ illustrates the temperature dependence of the 
decoherence rates for the case of one common bath. The 
values of the decoherence rates for the two bath case dif- 
fer insignificantly from the single bath case. We observe 
that the spread of the magnitude of the different deco- 
herence rates increases at intermediate mixing angles. As 
expected, the magnitude of the decoherence rates is max- 
imum in the case where the system Hamiltonian and the 
coupling to the bath are perpendicular to each other. For 
the opposite case, where the system Hamiltonian and the 
coupling to the bath commute, the decoherence rates van- 
ish for decreasing temperature, i.e., only flipless dephas- 
ing processes contribute to the overall decoherence. Note 
that in the case where the system Hamiltonian and the 
coupling to the bath commute (6\ — tt/2 and 62 = tt/2) 
the relaxation rate vanishes. 

It is found that the dephasing rates depend strongly on 
the qubit asymmetry. Nevertheless the parameter spread 
of the qubit energies can be quite large (around 10%) 
without affecting the decoherence properties consider- 
ably. However, for very large asymmetries and a bath 
coupling, which is perpendicular to the system Hamilto- 
nian, the decoherence rates increase exponentially with 
asymmetry. 



IV. CONCLUSIONS 



A system of two pseudospins coupled by an Ising type 
zz-interaction, which models e.g. superconducting charge 
or flux qubits, was investigated. It was shown that the 
for the system of two pseudospins the optimum gate per- 
formance of different gate operations is closely related to 
their composition of elementary gates and the coupling 
to the bath. When considering the gate quality factors, 
the temperature and special symmetries of the system- 
bath coupling have the largest influence on the decoher- 
ence properties, whereas the single or two bath scenarios 
make only little difference. 

For the CPHASE operation at low temperatures, the 
optimum gate quality factors are at pure a z system-bath 
coupling due to the fact that only in this case all indi- 
vidual Hamiltonians necessary for performing the gate 
operations and the system bath coupling commute. Sim- 
ilarily, the CNOT gate operation approaches the best gate 
quality factors close to a z system-bath coupling with a 
slight a x admixture due to the Hadamard operations. 

For very large temperatures, the temperature effec- 
tively symmetrizes the system and thus entanglement is 
always preserved during a CNOT gate operation indepen- 
dently of the system-bath coupling. 

Furthermore, it is found that the parameter spread of 
the tunnel matrix elements of the qubits, when operated 
close to the degeneracy point can be quite large (approx- 
imately 10%) for a bath coupling which commutes with 
the system Hamiltonian without affecting the decoher- 
ence properties considerably. In this case the differences 
in the decoherence rates stay below 1%. 
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